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Abstract 

We consider the change of entanglement of formation AE produced by a 
unitary transformation acting on a general (pure or mixed) state p describing 
a system of two qubits. We study numerically the probabilities of obtaining 
different values of AE, assuming that the initial state is randomly distributed 
in the space of all states according to the product measure recently introduced 
by Zyczkowski et al. [Phys. Rev. A 58 (1998) 883]. 
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Entanglement is one of the most fundamental issues of quantum theory [1]. It is a 
physical resource, like energy, associated with the peculiar non-classical correlations that 
are possible between separated quantum systems. Recourse to entanglement is required 
so as to implement quantum information processes [2-9] such as quantum cryptographic 
key distribution [10], quantum teleportation [11], superdense coding [12], and quantum 
computation [13-15]. Indeed, production of entanglement is a kind of elementary prerequisite 
for any quantum computation. Such a basic task is accomplished by unitary transformations 
U (quantum gates) representing quantum evolution acting on the space state of multipartite 
systems. U should describe nontrivial interactions among the degrees of freedom of its 
subsystems. How to construct an adequate set of quantum gates is one of the fundamental 
questions about quantum computation. A pleasing answer is found. Any generic two-qubits 
gate suffices for universal computation [16]. One would then be interested in ascertaining 
just how efficient distinct C/'s are as entanglers. Quite interesting work has recently been 
performed to this effect (see, for instance, [17-21]). 

A state of a composite quantum system is called "entangled" if it can not be represented 
as a mixture of factorizable pure states. Otherwise, the state is called separable. The above 
definition is physically meaningful because entangled states (unlike separable states) cannot 
be prepared locally by acting on each subsystem individually [22,23]. 

A physically motivated measure of entanglement is provided by the entanglement of 
formation E[p] [24]. This measure quantifies the resources needed to create a given entangled 
state p. That is, E[p] is equal to the asymptotic limit (for large n) of the quotient m/n, 
where m is the number of singlet states needed to create n copies of the state p when the 
optimum procedure based on local operations is employed. The entanglement of formation 
for two-qubits systems is given by Woofers' expression [25], 




(1) 



where 




(2) 
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and C stands for the concurrence of the two-qubits state p. The concurrence is given by 

C = max(0, Ai — A 2 — A3 — A4), (3) 

Aj, (i — 1, . . . 4) being the square roots, in decreasing order, of the eigenvalues of the matrix 
pp, with 

p = (a y ®a y )p*(ay®a y ). (4) 

The above expression has to be evaluated by recourse to the matrix elements of p computed 
with respect to the product basis. 

In the present effort we will concern ourselves with one of the basic constituents of any 
quantum processing device: quantum logic gates, i.e., unitary evolution operators U that 
act on the states of a certain number of qubits. If the number of such qubits is m, the 
quantum gate is represented by a 2 m x2 m matrix in the unitary group U(2 m ). These gates 
are reversible: one can reverse the action, thereby recovering an initial quantum state from 
a final one. We shall work here with m — 2. The simplest nontrivial 2-qubits operation 
is the quantum controlled-NOT, or CNOT (equivalently, the exclusive OR, or XOR). Its 
classical counterpart is a reversible logic gate operating on two bits: e±, the control bit, and 
e 2 , the target bit. If ei = 1, the value of e 2 is negated. Otherwise, it is left untouched. 
The quantum CNOT gate C 12 (the first subscript denotes the control bit, the second the 
target one) plays a paramount role in both experimental and theoretical efforts that revolve 
around the quantum computer concept. In a given orthonormal basis (|0), |1)}, and if we 
denote addition modulo 2 by the symbol ©, we have [26] 

|ei) |e 2 ) -> C12 -> |ei) |ei © e 2 ). (5) 

In conjunction with simple single-qubit operations, the CNOT gate constitutes a set of gates 
out of which any quantum gate may be built [16]. In other words, single qubit and CNOT 
gates are universal for quantum computation [16]. 

As stated, the CNOT gate operates on quantum states of two qubits and is represented 
by a 4x4-matrix. This matrix has a diagonal block form. The upper diagonal block is just 
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the identity 2x2 matrix. The lower diagonal 2x2 block is the representation of the one-qubit 
NOT gate U NO t, of the form 

1 

1 (6) 
A related operator is Ug, for which the lower diagonal block is of the form 

cos 9 sin 9 

- sin 9 cos 9 (7) 

C12 is able to transform factorizable pure states into entangled ones, i.e., 

C12 : [ci|0> + c 2 |l)]|0) <- Cl |0)|0) + c 2 |l)|l), (8) 

and this transformation can be reversed by applying the same CNOT operation once more 
[26]. In general, the action of Ucnot ( Ug) on a 2-qubits state (pure or mixed) produces a 
change of entanglement AE. 

The two-qubits systems with which we are going to be concerned here are the simplest 
quantum mechanical systems exhibiting the entanglement phenomenon and play a funda- 
mental role in quantum information theory. The concomitant space 7i of mixed states is 
15-dimensional and its properties are not of a trivial character. There are still features of 
this space, related to the phenomenon of entanglement that have not yet been character- 
ized in full detail. One such characterization problem will be the center of our attention 
here. We shall perform a systematic numerical survey of the action of Ucnot ( Ug) on 
our 15-dimensional space in order to ascertain the manner in which P(AE) is distributed 
in 7i, with P the probability of generating a change AE associated to the action of these 
operators. This kind of exploratory work is in line with recent efforts towards the systematic 
exploration of the space of arbitrary (pure or mixed) states of composite quantum systems 
[27-29] in order to determine the typical features exhibited by these states with regards to 
the phenomenon of quantum entanglement [27-33]. It is important to stress the fact that we 
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are exploring a space in which the majority of states are mixed. The exciting investigations 
reported in [17-21] address mainly pure states. 

As an illustration of the type of search we intend to perform, consider the real pure state 

p = a \00) + b\01) + c|10) + d\ll), (a, b, c, d E K + ), a 2 + b 2 + c 2 + d 2 = 1, (9) 

whose concurrence squared is 

C 2 = 4{ad - be) 2 , (10) 

and the transformations 

p' = UcNOT pUcNOT 
f> = U*/2 f> U~f 2 

P = U«/4pUjf (11) 

The ensuing squared concurrences for p' are, respectively, 

Ccnot — 4(ac — bd) 2 
Cl /2 = 4(ac + bd) 2 

C 2 /4 = 2(a 2 + b 2 )(c 2 + d 2 ) + 4[(-a 2 + b 2 )cd + (c 2 - d 2 )ab}. (12) 

Different changes of entanglement are seen to take place. The question is: given an initial 
degree of entanglement of formation E, what is the probability P(AE) of encountering a 
change in entanglement AE upon the action of Ucnot (Ug)7 

Due to the relevance of the Ucnot gate, we are going to focus our present considerations 
upon the probability distribution P(AE) associated with Ucnot- However, we must bear 
in mind that the P(AE)'s associated with other gates may be different from the P(AE) 
generated by the Ucnot gate. As an illustration of this differences we also studied some fea- 
tures of the P(A_E)'s corresponding to members of the mono-parametric family of gates Ug. 
This family of unitary operations comprises as a particular member the identity operation 
on two-qubits, / = Ug =0 , for which we trivially have P(AE) = for all AE ^ 0. 
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To answer the type of questions mentioned above we will perform a Monte Carlo explo- 
ration of Ti. To do this we need to define a proper measure on Ti. The space of all (pure and 
mixed) states p of a quantum system described by an N- dimensional Hilbert space can be 
regarded as a product space S = V x A [27,28]. Here V stands for the family of all complete 
sets of orthonormal projectors {Pi}f =1 , J2iPi — I {I being the identity matrix). A is the 
set of all real ./V-tuples {Ai, . . . , \n}, with < A« < 1, and J2i K — 1- The general state in 
S is of the form p = X^AjPj. The Haar measure on the group of unitary matrices U(N) 
induces a unique, uniform measure v on the set V [27,28,34]. On the other hand, since the 
simplex A is a subset of a (N — l)-dimensional hyperplane of TZ N , the standard normalized 
Lebesgue measure £at_i on TZ N ^ provides a natural measure for A. The aforementioned 
measures on V and A lead to a natural measure p on the set S of quantum states [27,28] , 

p = v x £jv_i. (13) 

We are going to consider the set of states of a two-qubits system. Consequently, our 
system will have N = 4 and, for such an N, S = TL. All our present considerations are 
based on the assumption that the uniform distribution of states of a two-qubit system is the 
one determined by the measure (13). Thus, in our numerical computations we are going to 
randomly generate states of a two-qubits system according to the measure (13) and study 
the entanglement evolution of these states upon the action of quantum logical gates Ug. 

Before embarking in our exploration of the entanglement changes produced by unitary 
operations, it is convenient to briefly review the salient features of the distribution of entan- 
glement on the state-space TC of two-qubits. Fig. la plots the probability P(E) of finding 
two-qubits states of Ti endowed with a given amount of entanglement of formation E. The 
solid line corresponds to all states (pure and mixed), while the dashed curve depicts pure 
state behavior only. We clearly see that our probabilities are of a quite different character 
when they refer to pure states than when they correspond to mixed ones. Most mixed states 
have null entanglement, or a rather small amount of it, while the entanglement of pure states 
is more uniformly distributed, (see the enlightening discussion in [27]). 
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Fig. lb plots the probability P(E F ) of generating via the CNOT gate a final state (pure 
or mixed) with entanglement Ep if the initial entanglement is zero (solid line). The mean 
amount of CNOT-generated entanglement from an initial state with E = is (Ep) = 0.0052. 
For the 7r/4-gate this figure equals 0.0023. For the sake of comparison we plot alongside the 
solid line a dashed curve that depicts the probability P{E) of randomly selecting a state 
endowed with an entanglement E. We encounter a mean entanglement (E) = 0.03. We 
appreciate the fact that it is quite unlikely that we may generate, via the CNOT gate, a 
significant amount of entanglement if the initial state is separable. 

Our gates act on an initial state of entanglement E and produces a final state of en- 
tanglement E F . In Fig. 2 we study (for all states) how the mean final entanglement (E F ) 
depends on the initial entanglement E . The final states are generated via i) the CNOT 
(solid line), ii) the U{9 = 7r/3) (dot-dashed line), iii) the U{9 = 7r/4) (dashed line), iv) the 
U{9 = 7r/6) (dotted line), and v) the U{9 — 0) (thin dashed line) gates. On average, these 
gates tend to disentagle the initial state. Here we have, of course, the trivial exception of 
the identity gate, 9 = 0, which does not produce any change of entanglement at all. As can 
be expected, the "disentangling" behaviour of the U{9) gate (when acting on mixed states) 
becomes less and less important as 9 approaches zero. As shown in the inset, the situation 
is different for pure states. In this case the CNOT gate (solid line) increases the mean 
entanglement up to ~ 0.5 for states with E lying in the interval [0, 0.5] and (E F ) becomes 
moderately smaller than E for states with large entanglement. The result is similar using 
the 7r/4 gate (dashed line) but, for small initial entanglement, {Ep) monotonically increases 
with E , reaching the value 0.738 for E — 1. The crossing between CNOT and 7r/4 results 
takes place at E Q = 0.53. The mean entanglement {E) = [29] for pure states is also 
drawn (thin "horizontal" line). 

Fig. 3a is a P(AE) vs. AE plot for pure states. The CNOT gate is compared with the 
unitary transformation (7) with 9 = it/ 4. In both cases there is a nitid peak at AE = 0. 
Thus, if the initial state has entanglement E, our survey indicates that the most probable 
circumstance is that the gate will leave it unchanged. One appreciates the fact that P{AE) 
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decreases steadily as AE grows. In the case of the 7r/4 gate the whole of the horizontal 
[—1, l]-interval is not accessible. 

It is important to bear in mind that, due to the invariance of the product measure un- 
der unitary transformations, if the initial states are uniformly distributed according to that 
measure, the probability distribution P(E F ) associated with the entanglement of the final 
states is the same as the distribution P(E ) characterizing the entanglement of the initial 
states. Consequently, it is instructive to compare the distribution P(AE) associated with a 
unitary gate with the distribution Pr(AE) obtained generating the final state randomly and 
independently from the initial one (that is, instead of using a unitary gate to generate the 
final state we generate it in random fashion). The dotted-dashed curve in Fig. 3a depicts 
the probability Pr(AE) of obtaining a difference AE between the amounts of entanglement 
of two randomly selected pure states. This last curve serves as a reference pattern for appre- 
ciating the net gate effect, i.e., the difference between the unitary gate-governed evolution 
and that of a random change-of-state process. This global picture changes dramatically if 
we consider mixed states as well (see Fig. 3b). For both gates considered here, there is still 
a peak at AE = 0, but of a much sharper nature. 

We are going to consider next the set of initial states that, as a result of our gate 
operation, suffer a given, fixed entanglement's change AE. To that effect we depict in 
Figs. 4a and 4b the behaviour of the mean initial entanglement (E ) as a function of the 
entanglement's increase AE. In other words, for each possible value of AE, we calculate the 
mean value of the entanglement of formation associated with all those initial states yielding 
the (same) change of entanglement AE upon the action of the CNOT or the ir/4 gates. Fig. 
4a corresponds to pure states and Fig. 4b to all states. In Figs. 4c (pure states) and 4d 
(all states) we plot (Eq) — (E ) 2 vs. AE. Pure states behave again in a drastically different 
manner as that of mixed states. We appreciate in Fig. 4a the fact that, for pure states, 
the CNOT exhibits the same qualitative behaviour as the ir/4 unitary transformation. The 
interval [—1,1] of AE for the 7r/4 transformation is not wholly accessible. We see that if the 
average amount of entanglement is large, then the two types of action we are considering here 
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will tend to disentangle the initial state. If mixed states enter the scene a totally different 
picture emerges. At the origin (AE = 0), an appreciable change of curvature is noticed. 
Confirming the results of Fig. 3b, we see that if the initial entanglement is zero, it tends to 
remain so. The dispersion is largest at the origin, for pure states (Fig. 4c). Including all 
states (Fig. 4d) introduces once again dramatic effects, the dispersion is negligible at the 
origin, grows rapidly in symmetric fashion, attains symmetric maxima and then decreases 
steadily towards zero for large A_E-values. The dispersion is much larger for two states 
connected through the 7r/4 operation than for two states connected via CNOT. 

Fig. 5 depicts for the CNOT gate (left) and for the 7r/4-gate (right) the probability 
distribution P(E F ) vs. E F for fixed initial entanglement E = 0.1, 0.2, 0.3, and 0.4, respec- 
tively. Notice that several crossings take place. No monotonous behavior can be therefore 
be associated to the E^-change. 

There are several information measures that are in common use for the investigation of 
quantum entanglement. The von Neumann measure 

S 1 = -Tr(plnp), (14) 

is important because of its relationship with the thermodynamic entropy. On the other 
hand, the so called participation ratio, 

W = (is) 

is particularly convenient for calculations and can be regarded as a measure of the degree 
of mixture of a given density matrix [17,27,30]. It varies from unity for pure states to N for 
totally mixed states (if p is represented by an iV x N matrix). It may be interpreted as the 
effective number of pure states that enter the mixture. If the participation ratio of p is high 
enough, then its partially transposed density matrix is positive, which for N = 4 amounts 
to separability [35,36]. 

The so-called g-entropies, which are functions of the quantity 

u;, = Tr(p«), (16) 
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provide one with a whole family of entropic measures. In the limit q — > 1 these measures 
incorporate (14) as a particular instance. On the other hand, when q = 2 they are simply 
related to the participation ratio (15). 

We revisit next how i) the mean initial entanglement (E ), and ii) the associated disper- 
sion (Eq) — (E ) 2 , characterizing those initial states yielding a given increase in entanglement, 
behave as a function of the aforementioned entanglement's change AE. We consider the 
action of the CNOT gate upon states having a given amount of "mixedness" (as measured 
by the participation ratio). We shall consider two such values, namely, R — 1.4 (left side 
of Fig. 6) and R = 2.2 (right side of Fig. 6). As it is expected, the maximum value of 
the mean initial entanglement {E ) decreases as the degree of mixture (R) increases. The 
corresponding graphs are depicted in Figs. 6a and 6b, which are to be compared to those 
of Figs. 4b and 4d, respectively. Quite significant differences are observed, specially in the 
case of dispersion, where we have now maxima instead of minima at the origin AE = 0. 

Two-qubits systems are the simplest quantum mechanical systems exhibiting the entan- 
glement phenomenon and play thereby a fundamental role in quantum information theory. 
The properties of its associated, 15-dimensional space Ti of mixed states are of a highly 
non-trivial character. As a consequence, the entanglement-related features of H are the 
subject of continuous interest as they have not been characterized in full detail yet. In 
this work we have investigated one of these characterization problems: the workings of logic 
gates, and in particular of the controlled-NOT gate, an operation that establishes or removes 
entanglement of formation E. 

To such an effect we performed a systematic numerical survey of the action of Ucnot 
(and of more general gates Ug) on our 15-dimensional space. The underlying goal was 
that of ascertaining the manner in which the probability P(AE) of generating a change of 
entanglement AE, associated to the action of these gates, is distributed in 7i. 

We have found that the statistical characteristics of our gates are quite different when 
they operate on mixed states vis-a-vis their effect on pure states. 7i is heavily populated 
with separable mixed states. We have shown that the probability of entangling these states 
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with our logic gates is rather low. Also, acting on an entangle mixed state is more likely 
that the gates will diminish the entanglement rather than augmenting it. 
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FIGURE CAPTIONS 



Fig. 1-a) Probability of finding two-qubits states with a given amount of entanglement of 
formation E. The solid line corresponds to all states (pure and mixed) and the dotted line to 
pure states, b) Probability of generating via the CNOT gate a final state with entanglement 
E F if the entanglement E of the initial state is zero (solid line). The dashed line is the 
probability P{E) of randomly selecting a state with entanglement E. 

Fig. 2- Mean final entanglement {Ep) vs. its initial value Eq for all states obtained via the 
CNOT gate (solid line), the U{9 = tt/3) (dot-dashed line), the U{9 = tt/4) (dashed line), the 
U(9 = 7r/6) (dotted line), and the U(9 — 0) (thin dashed line) gates. The inset illustrates 
the concomitant results for pure states. 

Fig. 3- a) P(AE) vs. AE for pure states. The change of entanglement AE arises as a result 
of the acting of a CNOT gate (solid line) and a 7r/4-one (dashed curve). The dotted-dashed 
curve reflects the entanglement change AE between two randomly chosen pure states, b) 
Corresponding results for all states. The solid line corresponds to the CNOT transformation 
and the dashed one to the 7r/4 gate. 

Fig. 4-a) Mean initial entanglement (E ) vs AE for pure states, b) For all states, c) 
{El)- {Eq) 2 vs. AE for pure states, d) For all states. The solid line refers to the CNOT 
transformation and the dashed one to the 7r/4 gate. 

Fig. 5. Effects of the action of the CNOT gate (left) and of the 7r/4-gate (right). We plot 
P(E F ) vs. E F for fixed initial entanglement E = 0.1, 0.2, 0.3, and 0.4, respectively. 

Fig. 6 a) Mean initial entanglement (Eq) vs AE via the CNOT gate for states with a given 
value of the participation rate R = 1.4 and R = 2.2. b) Its associated fluctuation {Eq)~ 
{E ) 2 vs. AE. 
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